Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN KTeX style file v2.2) 



Emission-line Diagnostics of Low Metallicity AGN 



Brent A. Groves \ Timothy M. Heckman^ and Guinevere Kautfmann^ 

^ Max-Planck Institut fur Astrophysik, Karl-Schwarzschild Strasse 1, Garching D-85748, Germany 



Dept. of Physics and Astronomy, John Hopkins University, Baltimore, MD 21218, USA 



o 
o 



5 February 2008 



o 
o 

o 



ABSTRACT 

Current emission-line based estimates of the metallicity of active galactic nuclei (AGN) at 
both high and low redshifts indicate that AGN have predominantly solar to supersolar metal- 
licities. This leads to the question: do low metallicity AGN exist? In this paper we use pho- 
toionization models to examine the effects of metallicity variations on the narrow emission 
lines from an AGN. We explore a variety of emission-line diagnostics that are useful for iden- 
tifying AGN with low metallicity gas. We find that line ratios involving [Nii] are the most 
robust metallicity indicators in galaxies where the primary source of ionization is from the 
active nucleus. Ratios involving [S ii] and [Oi] are strongly affected by uncertainties in mod- 
elling the density structure of the narrow line clouds. To test our diagnostics, we turn to an 
analysis of AGN in the Sloan Digital Sky Survey (SDSS). We find a clear trend in the relative 
strength of [N n] with the mass of the AGN host galaxy. The metallicity of the ISM is known 
to be correlated with stellar mass in star-forming galaxies; our results indicate that a similar 
trend exists for AGN. We also find that the best-fit models for typical Seyfert narrow line 
regions have supersolar abundances. Although there is a mass-dependent range of a factor of 
2-3 in the NLR metallicities of the AGN in our sample, AGN with sub-solar metallicities are 
very rare in the SDSS. Out of a sample of ~23000 Seyfert 2 galaxies we find only ~40 clear 
candidates for AGN with NLR abundances that are below solar 
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1 INTRODUCTION 



In recent years there has been increasing evidence that the growth 
of supermassive black holes at the centres of galaxies is closely 
linked to nuclear star formation and the formation of spheroids 
(eg ^ichstone et al. 1998; Haehnelt et al. 1998). The tight corre- 
lation between_the^ black hole mass and bulge velocity disper- 
sion jperrarese & Merrittii200Ck iGebhardt et al..200Q,') . as well as 
the w eaker correlation between black hole mass and sp heroid 
mass fKormendv & Richstondfl99S iMagorrian et al]ll998h . indi- 
cate that the fuelling of the central black hole must be accompanied 
by spheroid growth. Although part of the spheroid growth may re - 
sult from mergers of stellar disks (eg KaufFmann & Haehnelj2000l) . 
star formation in the central regions of the gal axy is also expected 
to contribute to the overall growth of the bulge iSchmitt et all 199^ 
Icid F ernandes et al. 2001, 2004). 

Active galactic nuclei (AGN) are found in galaxies where 
black holes are growing through the accretion of surrounding gas. 
The gas is also likely to be associated with the nuclear star forma- 
tion. By measuring the element abundances of the gas surrounding 
the AGN we obtain an indirect measurement of the star formation 
history of the host galaxy. If most of the gas has been converted 
into stars before the AGN becomes visible, then the gas is likely to 
be of high metallicity. If, on the other hand, the AGN is concurrent 



with or precedes the main episode of star formation in the bulge, 
low metallicity gas is likely to be found in some AGN. 

There have been analyses of elemental abu ndances in AGN 
for over two decades jPavidson & Netzeilll979h . In high-z AGN, 
this work has mainly focused on the broad emission lines to de- 
termine abundances. There have also been a number of studies 
that have used absorption line and narrow emission line diag- 
nostics (see Hamann & Ferland 1999, for a review). In nearby 
AGN the narrow emission lines are readily observable, and are 
oft en the ea siest met hod to determine m etallicities as pioneered 
by IStorchi B ergman n & Pastorizal il989h . The results for both 
high and low redshift AGN appear to show that most AGN 
have sola r to supersolar metallicities ( Storchi-Bergmann e t all 
[1998; Ha mann et ai]|2002h . Although the results do depend on 
the detailed model ass umptions (see eg iHamann et alj l2002t 
[Komossa & Schui3l993), there does appear to be some consensus 
that the gas in AGN is usually of high metallicity. 

The paucity of AGN with sub-solar metallicities may be linked 
to the fact that AGN in low mass galaxies appear to be very rare i n 
the local Universe, (see e.g.lGreene & Ho''2004':'Barth et al."2005l). 
Gas-phas e metallicity is known to be strongly correlated with stel- 
lar mass iTremonti et alj|2004l) . If AGN are predominantly found 
in massive, bulge-dominated galaxies that have processed most of 
their gas into stars by the present epoch, it may not be too surprising 
that the inferred metallicities are almost always high. 
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Although these conditions may hold at the present epoch, it 
is interesting to speculate whether the host galaxies of AGN might 
be considerably different at higher redshifts. For example, in the 
models of Kauffmann & Haehnelt (200G) the typical gas fractions 
in AGN hosts are expected to evolve from 5 - 10% at the present 
day, to ~ 20% at z = 1 and ~ 50% at z ~ 2. Even if high redshift 
AGN reside in massive galaxies, the gas-phase metallicities of their 
hosts might still be expected to be lower than at the present day. 

In this paper, we explore these questions by using photoion- 
ization models to examine the effects that a variation in metallic- 
ity would have upon the narrow emission lines from an AGN. We 
ask whether there are particular emission-line diagnostics that are 
most useful for identifying AGN with low metallicity gas. Using 
the models as a guide, we then search for low metallicity AGN 
in the Sloan Digital Sky Survey (SDSS). The SDSS provides high 
quality spectra from which we can measure both the host galaxy 
and the emission-line properties of a very large sample of nearby 
AGN, so it is an ideal database for searching for rare cases of low 
metallicity AGN in the nearby Universe. Finally, we discuss the 
implications of our results for searches for low metallicity AGN at 
higher redshifts. 
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Figure 1. The lBaldwin. PhilHDS. & TerlevichI <198lh emission line ratio di- 
agnostic diagram of [Nn]/i6584A/Hff versus [Oni]/i5007A/H/3 is plotted 
for SDSS emission line galaxies. Star-forming galaxies are shown in green, 
Seyfert 2s in blue and LINERs in red. The two diiferent curves used to 
separate AGN and star bursts are indicated by tlie solid curve (Kauifmann 
separator) and dashed curve (Kewley separator). 



2 NARROW LINE REGIONS AND LOW METALLICITY 

The use of nebular emission lines produced in Hii regions to in- 
fer the gas -phase abun dances of extragalactic objects has a long 
history (eg lAlleilfigil . Once understanding grew of the nuclear 
emission from AGN it was not long before the use of emission 
lines as abundance diagnosti cs was extend ed to these objects (eg 
iDavidson & NetzeJl979l:IStorchi Bergmann & Pastorizall989l) . 

There is a good understanding of the ionizing sources in H ii 
regions (the O <& B stars) and photoionization models to treat 
the physics of these regions are quite well-developed. As a re- 
sult, the correlation of emission line ratios with abundances in 
sta r-forming regions is reasona bly well understood. R ecent papers 
bv lCharlot & Longhettil J200lh and lKewlev & Dopital i2002>^ have 
provided calibrated relationships between gas-phase oxygen abun- 
dance and combinations of strong emission lines that are observ- 
able at optical wavelengths. 

For AGN the situation is markedly different. The ionizing 
source is not as well understood, and has a much greater range in 
luminosity. The ionizing spectrum is also harder than that produced 
by stars, and shocks are also more likely to contribute to the ion- 
ization state. While recent theoretical models (Baldwin et al. 1995^ 
iGroves etal1 E004a) have helped clarify some of the ionization 
structure, considerable uncertainties still exist. Many AGN show 
both broad and narrow line emission, and the ionization conditions 
in both the broad and narrow emission line regions show strong 
variations between different AGN. 

Both the broad and the narrow emission lines can be used to 
estimate gas phase abundances. The broad line region (BLR) sam- 
ples the gas closer to the nucleus. However, these emission lines 
are only seen in Type 1 AGN , and the strong emission lines used 
for abundance analysis are only found in the UV. As a result, abun- 
dance analysis is restricted to high re dshift AGN or to objects w ith 
UV observations from space (see e.g. lHamann & Ferlantf 1999). 

The narrow line regions (NLR) of AGN are more readily ob- 
servable. Their greater spatial extent means they are less likely to 
be obscured. The NLR also produces strong optical emission lines 
from several species. The physics of the NLR is also simpler than 
that of the BLR; the physical conditions are in fact similar to H ii 



regions, so that many of the metallicity-sensitive emission line ra- 
tios that are routinely appli ed to star forming galaxies a re also good 
diagnostics in the NLR. .Storchi-Bergmann et'aP il998ii) have cali- 
brated several of these emission line ratios for nearby AGN, us- 
ing H II region determined metallicities and photoionization mod- 
els. The similarity between the NLR and the Hit regions in a 
galaxy also leads to ambiguities when estimating abundances. Un- 
less spatially-resolved observations can be obtained, it is often dif- 
ficult to disentangle the contributions to the emission line spectrum 
from the NLR and from star-forming regions within the galaxy. 

Figure [H shows one of the standard line ratio diagrams 
iBaldwiii! Phillips. & TerlevichI Il981l) used to distinguish star- 
forming galaxies, Seyfert 2 galaxies (or AGN NLRs) and Low Ion- 
ization Nuclear Emission-line Region galaxies (LINERs). A de- 
tailed description of the location of these different classes of ob- 
ject in the different line ratio diagrams and the criteria that de- 
fine our classificatio n system has been given in a recent paper by 
lKewlevetal1j200d) . 

In this particular diagram, the emission line galaxies are dis- 
tributed in V-shaped morphology, with the star-forming galaxies 
located in the left branch and the AGN in the right branch. The 
AGN branch actually consi sts of two classes of object, the Seyfert 
2 galaxies and the LINERs jKewlev et alJ20 0dh. 

The distribution of the starforming galaxies along the hori- 
zontal [Nii]/16584A/Hq' axis is strongly correlated with the metal- 
licity within the starforming region and the mass of the galaxy; 
both increase with increasing [Nn]/Ha. The vertical axis of 
[Oiii]/i5007A/Hy8 is associated with the average ionization state 
and temperature of the photoionized gas in the emission line galaxy. 
AGN are offset from the star-forming galaxies on this axis because 
of their much harder ionizing spectrum. Those objects which lie at 
the extremities of the AGN branch are "pure" AGN. As the contri- 
bution of star-formation to the emission line spectrum increases, 
the strength of [Oiii]/Hy8 decreases. This AGN-starburst mixing 
sequence culminates at the fulcrum of the V, where starforming 
galaxies and AGN dominated by starformation become indistin- 
guishable. When we model NLR spectra, we consider only "pure" 
AGN. 
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Table 1. Solar abundance & metallicity scaling 



Element 



Abundance Depletion 



H 0.000 0.00 

He -0.987 0.00 

C -3.61 -0.30 

N -4.20 -0.30 

O -3.31 -0.24 

Ne -3.92 0.00 

Na -5.68 -0.60 

Mg -4.42 -1.15 

Al -5.51 -1.44 

Si -4.49 -0.89 

S -4.80 -0.34 

CI -6.72 -0.30 

Ar. -5.60 0.00 

Ca -5.65 -2.52 

Fe -4.54 -1.37 

Ni -5.75 -1.40 



3 LOW METALLICITY NLR SIMULATIONS 

We have run a series of narrow line region (NLR) models with 
a wide range in both metallicity and ionization parameter'. The 
models are based upon the dusty, radiation pressure dominated 
models of Grov es et a l. 1 2004a b), but use a more realistic ioniz- 
ing spectrum iGroves et alj|2006t) . While it is unknown whether 
low metallicity AGN will contain dust, the observations of dust in 
high-z objects indica tes that such an assumption is not unreasonable 
fSertoldi et al. 2003). Unlike previous work on NLR metallicity (eg 
Istor chi-Bergmann et al, 1998), these models assume a more phys- 
ical structure for the NLR clouds and concentrate upon the effects 
of low metallicity. 

As the models have been discussed in depth in previous pa- 
pers, we only describe the relevant parameters and the changes that 
were made in order to address the issue of low metallicity AGN. 
Throughout this work we assume a fiducial value of «h ~ 10^ cm"' 
for the density as found in previous papers. 



3.1 Metallicity and Dust Depletion 

T he assumed abundanc es are based upon the solar abundance set 
of lAsplundetalJj2005h . which takes into account recent revisions 
of the solar abundances of several important elements like carbon 
and oxygen. Note that abundance ratios may vary between different 
AGN and even within the NLR of a single AGN; the abundances 
given in tableQare t hus meant to b e representative. 

As discussed in lGroves et alj j5)04a), the abundance with re- 
spect to hydrogen of most of the elements listed in table Q will 
scale linearly with total metallicity. The exceptions are He and N. 
For helium, the yield from stars must be added to the primordial 
abundance. In thi s work, we use the primordial measurements of 
[Pasel etan <1992h . which gives 



He/H = 0.0737 + 0.0293Z/Zo. 



(1) 



' The ionization parameter is a measure of the density of ionizing photons 
relative to the gas density. 

^ All abundances are logarithmic with respect to Hydrogen 
3 Depletion given as log(X/H)gas - log(X/H)isM 



As nitrogen possesses both primary and secondary nucleosyn- 
thesis components, the relationship between nitrogen abundance 
and metallicity differs from that of other elements. For low metal- 
licity galaxies (log (O/H) < -4.0) the N/0 ratio is approximately 
constant, as expected for a "primary" element whose production is 
independent of metallicity. However, for higher metallicity galax- 
ies (log (O/H) > -3.5) the N/0 ratio is found to rise steeply with 
metallicity. This suggests that nitrogen becomes dominated by sec- 
ondary producti on from CN O nucleosynthesis. 

We follow loroves et a l. (2004a) and use a linear combina- 
tion of the primary and secondary components of Nitrogen. This 
relatio nship is fitted to the data fromJV Iouhcine & Contini ( 200^) 
and Kennic utt. Bresolin. & Garnetl|j2003l) . with the requirement of 
matching the solar abundance patterns (see figure 2 in .Groves et all 
^004a). With the new abundances in tableQwe obtain the following 
relation. 



(N/H) = (O/H)(l0"'* + io(2"+tof=iolo/H])^ 



(2) 



Using these prescriptions, we have then explored six NLR 
metallicities: 4Zo, 2Zo, IZq, O.SZq, 0.25Zo, 0. IZq, and O.OSZq. This 
range is wide enough to explore trends in emission line proper- 
ties associated with a decrease in AGN metallicity. In terms of 
12 + log(0/H), these abundances correspond to 9.29, 8.99, 8.69, 
8.39, 8.09, 7 69 and 7.39 respectively. Note that with the new 
lAsplund et alJ (2005) abundances, the oxygen abundance at solar 
is approximately 0.2 dex lower than that given in previous work 
( Grevesse & Sauval 1998). As some of the metals present in the 
NLR clouds are depleted onto dust, the actual metallicity of the gas 
component in our models is approximately half of the total metal- 
licity. 

3.2 Narrow Line Region Dust 

While the total dust to gas ratio is known to decrease with de- 
creasing metallicity, the actual variation of dust abundance with 
metallicity is uncertain. In AGN no direct measurements can be 
made of dust depletion so for simplicity, we assume a solar de- 
pletion pattern for the dust and one that is constant with metal- 
licity. In Galactic regions with different physical conditions, the 
dust to gas ratio, and hence the depletion factors, are observed to 
vary iSavage & Sembachll996h . However, when regions of differ- 
ing metallicities are observed, the depletion factors appear to be 
approximately constant ( Vladilo 2002). Our use of constant deple- 
tion factors for NLR models of differing metallicities is therefore 
not unreasonable. 

The de pletion factors, given in table Q are based upon the re- 
cent work of lKimura et alj i2003{) who examined the metal absorp- 
tion along several lines of sight within the local bubble. The gas 
abundances were compared with solar abundances to determine the 
depletion. Some elements were not discussed in this work and fo r 
these we use previously assumed values f romlPopita et alj i200(l) . 
For all other dust properties we follow the'Groves et al.'feOO^ de- 
scription of NLR dust. Note that any variation of the depletion fac- 
tors will have a similar effect on the model to abundance variations. 



3.3 Ionizing Spectrum 

We follow Groves et alj i200d) and use an empirical fit to the 
lElvis et ^ ^gjobservations. 



/„ = v«Huvexp(-^)exp(-«iM) 
+ flv-xexp(-4)exp(-^^). 



(3) 
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Table 2. Input parameters for the models. The total incident flux density /q 
is scaled by 38.098 erg cm ^s-l, while Pogas/*: is scaled by lO^K cm . 



model# 


h 




So 


log(f/o) 


1 


4.0 


1.00 


20.16 


-0.40 


2 


3.0 


3.20 


4.73 


-1.03 


3 


2.0 


5.40 


1.87 


-1.43 


4 


1.0 


7.60 


0.66 


-1.88 


5 


0.5 


8.70 


0.29 


-2.24 


6 


0.25 


9.25 


0.14 


-2.57 


7 


0.1 


9.58 


0.05 


-2.98 



with ttEuv = -1.75, kTuv = 120 eV, kTsBB = 7.0 eV. 

The X-ray part of the spectrum has an index of ax = -0.85 
and an upper cut-off of kTx = 10' eV. The parameter a is set to 
a = 0.0055, which results in an optical-X-ray slope of index q-q-x ~ 
-1.4. For simplicity, we use the same ionizing spectrum for AGN 
of all metallicities, but we note that a lower abundance of metals 
and especially dust is likely to alter the structure of the AGN torus 
and possibly even the accretion disk. 

3.4 NLR model Parameters 

For each metallicity we mode l a range of NLR i onization condi- 
tions. We follow the work of Groves et alj i2006h and maintain a 
constant total pressure across the models, while varying both the 
incident flux density and initial gas pressure. The total pressure we 
consider is Ptot/k = lO^'Kcm"''. This corresponds to a [Sii] den- 
sity of rill - 10^ cm"^ which is a reasonable estimate of the NLR 
number density. While the effect of radiation pressure upon dust 
becomes less important at low metallicities, we maintain the same 
input parameters across the models for comparison purposes. 

In table|2|we describe the parameters used for the model sets, 
with the same inputs used for each metallicity. For each input flux 
density and initial gas pressure we give the resulting initial i on- 
izatio n state in terms of two parameters. Ho and Do ( Groves e t alj 
l2006h . Ho is a measure of the radiation pressure versus gas pressure 
(Hq = lolicPgas)), while Uq measures the relative number of ion- 
izing photons 5., to a normalized initial density {Of, = S,l{chii), 
fin = Pg-Jikmi 

The parameters explored here span the range in ionization 
conditions that are observed in typical narrow line region clouds. 
The models are all truncated at a column densityjrf log(A'(H i)) = 
2 1 , which is reasonable for NLR clouds fe.g. jCrenshaw et aL2003.") . 



4 METALLICITY EFFECTS ON THE EMISSION LINES 

When the total metallicity of a photoionized nebula decreases, there 
are several effects that cause changes in the final emission line 
spe ctrum. Many of these effects have b een described elsewhere 
(eglOsterbrockll989HDopita & Sutherlanci2003l:iKewlev & Dopital 

\2om . 

The dominant effect is the decreased abundance of an emitting 
species relative to hydrogen. This leads to weaker emission in lines 
such as nitrogen and oxygen relative to hydrogen recombination 
lines such as Ha and Hj8. The line ratios relative to the other abun- 
dant primordial element, helium, also change in a similar way. This 
effect is not as strong, because helium also decreases with metal- 
licity. 



1,0 



0.5 - 



-0.5 



-1.0 - 



L0.1Z. 2"; 



2Z. 



1Za 



0.5Z. 
0.26Z. i , . 




-2.5 



-2,0 



-1.5 



-1,0 
log([Nll]/Ha) 



-0.5 



0,0 



0.5 



Figure 2. The BPT diagram of [N n]/Ha versus [O in]/HjS demonstrating 
the decrease of metal emission lines relative to hydrogen. This decrease 
is stronger for the nitrogen lines due to the secondary component in this 
element. Model metallicities are as labelled, with the small numbers next 
to model indicating model run and associated ionization parameter as given 
in table |2] The background image shows the log scale density distribution 
of the SDSS emission line sample, with Seyfert galaxies lying in the right, 
blue colored branch. Each branch has been scaled separately to emphasize 
the position of the Seyfert galaxies. 



At metallicities above 0. IZg, nitrogen is dominated by its 
secondary component. A decrease in metallicity down to this 
value results in a decrease of the strength of the nitrogen emis- 
sion line relative to that of other metals. Below Z ~ O.IZq, 
the primary component is dominant and the nitrogen tracks the 
behaviour of the other metal lines more closely. These trends 
are cle arly seen in figure ID which p resents the models on the 
classic Baldwin . Phillips. & Terl evich ('1981") (BPT) diagram of 
[Nn]i6584A/HQ- versus [O in]/15007A/Hye. This diagram clearly 
reveals the strong decrease of the nitrogen emission lines relative 
to hydrogen as the metallicity decreases. 

As metallicity decreases, the temperature of the nebula also 
increases. This is a result of the decrease in the efficiently cool- 
ing metal emission lines and the availability of more hi gh en- 
ergy photons to ionize hydrogen (see Sutherland & Dopital [l993l 
for detailed explanations). Higher nebular temperatures affect 
the strength of temperature sensitive emission lines, such as 
[O m]/i5007A, which increases with temperature. Line ratios like 
[O m]/i5007A/Hj8 will thus decrease at a relatively slower rate with 
metallicity. This is also clear from the model presented in Figure 2. 
While temperature sensitive line ratios change with the metallicity 
of the gas, they are also strongly affected by the ionization state of 
the gas. This means that shocks and the varying contribution of star 
formation and AGN to the ionization become important. 

Another effect of lowering the metallicity is that the opacity 
of the gas is reduced, meaning a greater volume is needed to ab- 
sorb all ionizing photons. In AGN this effect causes noticeable dif- 
ferences because of the stronger contribution to the ionizing flux 
from X-rays, for which the hydrogen opacity is small. These X-rays 
give rise to the partially ionized region in NLRs, where important 
optical emission lines such as [Oi]6300A and [S n]6713-H31A are 
emitted. Models with lower metallicity have larger partially ion- 
ized re gions relative to the fully ionized region (see lGroves et alJ 
l2004al) . In our column density limited models, the partially ionized 
region becomes truncated at low metallicities. As a result, we are 
not confident in our ability to use these models to calibrate some 
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of the metallicity sensitive ratios commonly used in determining 
the abundances in star-forming galaxies, such as [Nii]/[Sii] and 
[S n]/Ht!' for use in the determination of AGN abundances. 



5 DIAGNOSTIC DIAGRAMS FOR LOW METALLICITY 
AGN 

Generally, individual AGN are best reproduced with a combination 
of models, varying not only the metallicity, but also the ionizing 
spectrum, the density, the incident flux or ionization parameter and 
even the geometry (see eg Oliva et al. 1999). 

However when looking at large numbers of AGN the best 
method is to look for trends or possible relation ships between 
the line ratios and metallicity (eg I Storchi-Bergma rm et aljfl998l : 
iNasao et al."2Q02': ' Groves et al .' 2004b). Previous AGN modelling 
has found that AGN as a group are best fit with models of around 
2 - 4Zo, or 12 + log(0/H) ~ 9.0 - 9.3. Thus any AGN NLR found 
to have metallicities below solar could possibly be considered "low 
metallicity". 

As discussed in the previous section, line ratios relative to hy- 
drogen are good diagnostics for low metallicity AGN. This is es- 
pecially true for nitrogen with its secondary component as shown 
in figure |5| The diagnostic of [Nii]/16584A/Ha has been used as 
one of the main metallicity ind icators in AGN due to its sensitivity 
and the strength of th e lines ( S torchi Bergmann & Pastorizall989t 
IStorchi-Bergmannet aL.1998,) . In general, low metallicity AGN 
will lie to the left of the Seyfert branch on the Nn BPT diagram. 
The decrease in the [O iii]/Hj8 ratio is not as strong because of the 
temperature sensitivity of this ratio. 

Another possible diagram for metallicity diagnosis is the 
[Nn]^6584A/[0 n]^3727A versus [O m]i5007A/ [O ii].l3727A di- 
agram, shown in figure This again relies upon nitrogen to sep- 
arate the different metallicities, but also uses the temperature sen- 
sitive ratio [O iii]/[0 ii] as a secondary diagnostic. The [N ii]/[0 n] 
ratio will separate AGN with lower N ii relative to O ii. As N ii and 
O II are both low ionization species , the ratio will only be weakly 
alTected by the ionization or density conditions. In addition, star 
bursts and AGN have similar [N ii]/[0 ii] ratios for similar metallic- 
ities, so the abundance determination of an active galaxy is only 
slightly perturbed by any contribution from star formation. The 
main disadvantage of this ratio is that [Nii] and [On] are widely 
separated in wavelength and the ratio is hence much more sensi- 
tive to the effects of reddening. The [O iii]/[0 ii] ratio is similar to 
[Oiii]/H/? as a temperature sensitive diagnostic in that is good in 
separating the Seyferts from the star- forming galaxies. It also has 
the additional benefit of not being strongly sensitive to metallicity. 

Both these diagnostic diagrams rely heavily on the [N ii] line 
to estimate the metallicity. This means that the diagnostic is sensi- 
tive to abundance variations relative to solar within the NLR due to 
difi^ering star formation histories. 

At higher redshift, the [Nii] line is no longer observable, but 
UV di agnostics become useful (Nagaoetal. 2006; Groves et al. 
l2004d) . One possible diagram is shown in figure |31 which plots 
the Near-UV line ratios of [Oii]/i3727, 9A /[Neiii]i3869A versus 
[Nev]/i3426A/[Neiii]/t3869A. This diagram relies upon the tem- 
perature sensitivity of the [Ne iii] line to distinguish different metal- 
licities in the [O ii]/[Ne iii] ratio. The [Ne v]/[Ne iii] ratio, while 
weakly dependent upon metallicity, is very sensitive to the ionizing 
spectrum and ionization state of the gas. Thus this ratio provides a 
way of distinguishing star-forming galaxies, which have very little 
[Nev] emission, from AGN. As in figure |3| the strength of figure 
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Figure 3. Metallicity diagnostic diagram of [N ii]/t6584A/[0 ii]/l3727A ver- 
sus [O iii]/i5007A/ [O ii]/i3727A with models increasing in metallicity from 
left to right (as labelled). The background image shows the distribution of 
the SDSS emission line galaxies, with Seyfert galaxies, colored blue, lying 
in the vertical branch and starbursts in the horizontal branch. Each branch 
has been scaled separately to emphasize the position of the Seyfert galaxies. 
The contribution by star formation to the AGN decreases the [O in]/[0 n] ra- 
tio. 
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Figure 4. Near-UV metallicity diagnostic diagram of 
[Oii]i3727,9A/[Neiii]i3869A versus [Ne v]i3426A/[Ne iii]i3869A. Pure 
star forming galaxies would lie below these curves at low [Ne v]/[Ne in] 
values. 

lilies in its use of strong lines and its ability to distinguish star 
forming galaxies from AGN. 



6 AGN SELECTION WITHIN SDSS 

Within this work we use the SDSS Data Release 4 (DR 4) spec- 
troscopic galaxy sample <Adelman-McCarthv et ai]l200d) . which 
includes u-, g-, r-, i- and z-band photometry and spectroscopy of 
over 500,000 objects. The spectra are taken using 3-arcsec diame- 
ter fibres, positioned as close as possible to the centres of the tar- 
get galaxies. The flux- and wavelength-calibrated spectra cover the 
range from 3800 to 9200A, with a resolution of R -1800. 

At the median redshift of the sample (z ~ 0.1) the spectro- 
scopic fibre typically contains 20 to 40 percent of the total galaxy 
light, thus contain a component due to the host galaxy as well as 
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the AGN. As described in lTremonti et alJ i2004b . we subtract the 
contribution of the stellar continuum from each spectrum, using 
the best fitting combination of template spectra from the population 
synthesis code of Bruzual & Chariot ( 2003). The best-fitting model 
also places c onstraints on the star formation history and metallicity 
of the galaxv iGallazzi et alj2005l) . and can be used to estimate stel- 
lar masses and star-formation histories jKaufFmann et al. 2003a). 

W e have used the same sample criteria as in Kewlev et al. 
i2006h to extract our narrow-line AGN subsample. The sample is 
limited to redshifts above 0.02. We also apply a signal-to-noise cut 
{S IN > 3) on the six dominant emission lines; H/3 /14861A, [Om] 
ASOOTA, [Oi] i6300A, Ho- ^6563 A, [Nii] i6584A, and the dou- 
blet [Sii] /i6716, 31A. This gives a s ample of ~ 170,000 emission 
line galaxies. The lKauffmann et alj i20033) Nii empirical relation 
is then used to define our AGN sample. 



log([OIII]/Hj8) > 0.61/(NII]/Hq' - 0.05)) + 1.3. 



(4) 



This results in a sub-sa mple of ~ 50000 galaxies c lassified as AGN. 
The Nil dividing line of lKewlev & Dopitiji20o3) . 



log([OIII]/Hj8) > 0.61/(NII]/Hci' - 0.47)) + 1.19, 



(5) 



is used to select a subsample of 20,000 "pure" AGN where the ion- 
ization is dominated by the active nucleus rather than H ii regions. 

The AGN sample includes two classes of emission line ob- 
jects; Seyfert 2s and Low Ionization Nuclear Emission-line Re- 
gions (LINERs). In this work we consider only Seyfert galaxies, 
as our diagnostic models for the NLR only apply to these objects. 
To distinguish LIN ERs from Seyferts, w e use the recent empirical 
dividing lines from lKewlev et all i2006l) 



log([OIII]/H8) > 1.361og([OI]/Ha) + 1.4, 



(6) 



Our final sample consists of 8800 pure Seyfert galaxies rang- 
ing in redshift from z ~ 0.02 to ~ 0.36, with a median of z ~ 0.1. 



6.1 Accounting for Reddening in the SDSS Sample 

Line ratios such as [Nii]/[On] have a large separation in wave- 
length and thus need to be corrected for reddening. We use the 
Hfy/H/? ratio and we assume dust-free values of 2.86 for star- 
forming galaxies and 3.1 for Seyferts. The majority of AGN in our 
sample will have ongoing starformation and their dust-free values 
of this ratio should lie betwe en these values. We th en assume the 
power- law slope of A-°^ froml Charlot & Falll iioool) to correct for 
the attenuation in each galaxy. 



7 LOW METALLICITY AGN IN SDSS 

Measurem ents of the average gas metallicities in starforming galax- 
ies (fig. 6. lTremonti etai]l2004h show that there is a strong trend 
of increasing metallicity with increasing stellar mass. This trend 
is also sgen in the average stellar metallicities in normal galaxies 
(fig. S. iGallazzi et all2005h . 

Figure |5| shows the mass distribution of the full emission- 
line galaxy sample and the AGN sample. AGN a re clearly biased 
towards higher masses teauflFmann et alJl2003d) and they range 
from ~ 10'" to ~ IO"^Mo. If AGN have similar metallicities to 
star-forming galaxies, they shou ld have abundances i n the range 
12 -I- log(0/H) ~ 8.95 to ~ 9.13 <Tremonti et al.l2004 eqn. 3), or 
metallicities of 2-3Zq. As shown in figure|5|this is also the metal- 
licity range for which our AGN photoionization models best repro- 




Figure 5. The Stellar Mass distiibution of the full emission line galaxy 
sample (solid line) and the AGN (including LINER) dominated galaxies 
{dashed line) in Solar masses. The peaks have been scaled by 11013 and 
2286 for the full and AGN sample respectively. 




Figure 6. Variation of the median angle of the Seyfert branch, (pNii, with 
the host stellar mass, M» in solai' masses. 0imii is the angle relative to the 
X-axis on the BPT [Nii] diagram using the empirical point [-0.5,-0.6] as 
the vertex (see figure 0. Only AGN-dominated Seyfert classified galaxies 
are considered. 



duce the emission line ratios of the typical Seyfert galaxies within 
our sample. 

Figure|2|shows that as the metallicity decreases, AGN are ex- 
pected to migrate to the left on the [O iii]/H/3 versus [N iil/Ha BPT 
diagram. Following . Kewlev et al.. l2006) . we choose an empirical 
base point on the [N n] BPT diagram ([-0.5,-0.6], as marked on fig- 
ure0 and use this as the vertex to determine the angle each AGN in 
our sample makes with the j:-axis. We label this angle ^nh. Figure 
|6|shows the variation of the median value of this angle as a function 
of the stellar mass of the host. This figure clearly reveals a corre- 
lation between the host stellar mass and metallicity. The change in 
(pNii from logM, ~ 9.5 to logM, ~ 11.5 corresponds to a change 
of about 0.2 - 0.3 dex in log(0/H) in the AGN models. 

We then use a cut in stellar mass to separate out a sample of 
low metallicity AGN for closer inspection. We select AGN hosts 
with log(M./Mo) < 10.0. This limit selects the tail end of the AGN 
mass distribution, while maintaining sufficient numbers for analy- 
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Figure 7. The BPT [Nii]/Hq- versus [Oiii]/H/3 diagram showing the dis- 
tribution of the low mass AGN (crosses). The background image shows 
the log scale density distribution of the full SDSS emission line sample as 
in figure I2] The solid and dashed curves show the Kauffmann and Kew- 
ley classification lines respectively. The boxed crosses highlight the lowest 
metallicity low mass Seyfe rt galaxies (se e figurelsl. The solid squares and 
circles represent the Shapley et alj 12005) and Erb et al. l2006 ) high red- 
shift star-forming galaxies respectively, and show their position relative to 
the low metallicity AGN. (see 
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Figure 8. The metallicity sensitive [Nn]/t6584A/[On]i3727A versus 
[O iii]/i5007A/[O ii]/13727A diagram showing the distribution of the low 
mass Kewley-selected AGN galaxies (crosses). The background image 
shows the logscale density distribution of the full SDSS emission line sam- 
ple. The Seyfert galaxies, marked as blue on the image, lie on the vertical 
branch at high [N n]/[0 n], while the star-forming galaxies lie on the mainly 
horizontal curve. LINERs exist at the dark clump on the star-forming branch 
at high [N n]/[0 n] . The dotted lines indicate our arbitrary metal & AGN 
cuts, with our final low mass, low metallicity AGN highlighted by squares. 



sis. This yields 1800 low-mass AGN, of which 650 lie above the 
Kewley line. Figure0shows the distribution of the low-mass AGN 
(both Seyferts and LINERs) on the standard BPT [N 11] diagram. As 
can be seen, the low mass objects, while scattered, tend to lie to- 
wards the left of the AGN branch. Many of them are strongly clus- 
tered close to the Kauffmann dividing line. To reduce the effects of 
confusion and contamination by star-formation, we concentrate on 
the Seyfert galaxies that lie above the Kewley line. This reduces the 
sample to ~ 350 galaxies. 

To estimate the metallicity of the low mass AGN we 
plot these objects in the [Nii]/16584A/[Oii]/13727A versus 
[O m]/i5007A/[O n]/l3727A diagram described in the previous sec- 
tion. This is shown in figure |S| The use of the Kewley cut means 
that the majority of starburst dominated AGN have been removed 
from the sample as can be seen from the scarcity of objects in the 
heavily shaded greyscale region. However there are still some ob- 
jects which appear to lie in the starburst branch. These galaxies may 
either have incorrectly measured line fluxes or incorrect reddening 
corrections (as is certainly the case for the lowest [O in]/[0 11] ob- 
ject). Once again the low mass Seyferts tend to lie to the left, low 
metallicity side of the main Seyfert branch. 

Figure 8 demonstrates is that even at the low mass end of the 
AGN distribution, the AGN gas metallicity is still around solar (ac- 
counting for dust depletion). We note that this is similar to the 
gas-phase me tallicities of star-form ing galaxies of the same mass 
(equation 3 in lTremonti et all2004h . 

To select the most extreme low metallicity AGN in our sample, 
we apply two cuts to exclude those objects that lie within the main 
AGN or starforming branches. These cuts are shown by the dotted 
lines in figure|S|and are: 

log([OIII]/[OII]) > 2.1 log([NII]/[OII]) + 0.85 
log([OIII]/[OII]) > -1.21og([NII]/[OII])-0.7. ^ ' 

These cuts select 50 of the 559 low mass AGN as our low- 
est metallicity candidates, and they are plotted as boxed crosses on 
figures and [8| Note that one object has multiple spectra, which 



leaves 47 unique candidates. Six of these objects lie at redshifts 
greater than 0.1, with at least one being a strongly starbursting 
galaxy. These low "metallicity" AGN lie clearly to the left of the 
AGN branch in the BPT [Nii] diagram (see figure^, and similarly 
lie on the left part of the AGN branch in both the BPT [S 11] and 
[O i] diagrams. The few low mass AGN that lie to the left of the 
AGN branch on figure^and are not selected by our cuts are actu- 
ally at too low a redshift (z S 0.03) for accurate measurement of 
the [0 11] line flux. 
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Figure 9. SDSS fibre spectra of six of the low mass, low metalUcity selected AGN. a) shows the spectrum of a starburst galaxy that falls in our selection criteria 
while b) - f) shows five representative low metallicity AGN. The difference between the starburst and AGN can be seen, as well as other AGN signatures, such 
as the strong [Hen]/i4686A line (unmarked line to the left of HjS), or old stellar continuum. 
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Figure 10. SDSS images of the six sample galaxies from figurel^with a) an 
example starburst galaxy and b) - f) low mass AGN. Each image shows the 
same area, with the bar in each frame indicating 5". The image centre, as 
marked by the lines, is the position of the 3" SDSS fibre centre. Note that 
for most AGN the 3" fibre encompasses a significant fraction of the host 
galaxy. 



As a final check we examined the spectra of these objects and 
a representative selection of these is shown in figure|9| One of these 
objects (figure the galaxy with the lowest [Nii]/[Oii] ratio, is 
a clear starforming galaxy with an apparent incorrectly estimated 
[Om]5007 strength. However this is the only clear pure starburst 
galaxy among the 41 objects. Some of the spectra of the other 40 
AGN candidates exhibit clear signatures of young stars as well as 
He II /i4686A emission. Others are dominated by older stars, and 
there are also a number of AGN with strong B aimer absorption 
lines. All objects have very weak [Nii]6584 and [Oi]6300 (ex- 
pected because of their low metallicities), but the strong [O iii]5007 
emission line and relatively strong [Sii] 6716,31 doublet classify 
these objects as AGN. One other point to note is the presence of 
the [O iii]/14363A line in many of the spectra, indicative of the high 
gas temperatures in these objects, and also suggestive of low metal- 
licities. 

As a complement to the SDSS fibre spectra, we also show im- 
ages of the same sample in figure fTol The starforming galaxy has 
a distinctly different morphology compared to the AGN, appearing 
diffuse and blue in color. The images of the low mass AGN indicate 
that their host galaxies are faint and relatively compact. They tend 





Figure 11. SDSS images of six random NLR dominated Seyfert galaxies. 
These galaxies have redshifts z = 0.05 ± 0.01, and all fall within the main 
Seyfert branch on line diagnostic diagrams. Markings the same as in figure 

m 

to be red, with a disky or lenticular structure surrounding a rela- 
tively bright nucleus. For comparison, in figure ll ll we show Seyfert 
galaxies that fall in the main branch on line diagnostic diagrams and 
have a median redshift similar to the low mass sample of z ~ 0.05. 
The most obvious difference is their size. The spiral structure in the 
host galaxy is also much more evident. However in some respects 
the low mass and high mass AGN are similar; both have a bright 
red nucleus/bulge, and diffuse, blue outer disk. 

7.1 The Lowest Metallicity AGN 

In the previous section, we searched for low metallicity AGN by 
looking for objects offset from the main AGN branch. However the 
diagnostic diagrams of section|5]indicate that the lowest metallicity 
objects could actually lie within or even below the star-forming 
branch. 

To explore this possibility, we selected from the SDSS sam- 
ple all objects with a signal-to-noise greater than 3 in Ha, H/J and 
[O iii]/i5007A, and looked at the metallicity sensitive line ratios 
such as [S ii]/Hq' [Ne iii]/H/? or [N ii]/[0 n] . The spectra, images and 
galaxy parameters of all objects with extremely low values of these 
ratios were examined. 

All of these objects appeared to be indistinguishable from low 
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mass, low metallicity starburst galaxies. They all had blue colours 
and high equivalent width emission lines. None showed obvious 
signs of AGN activity, such as the Hen /i4686A emission line or 
[Ne v]/i3426A emission in the higher redshift objects. We found no 
clear evidence for extremely low metallicity AGN. If such objects 
exist, they are either very rare they are indistinguishable from a low 
metallicity starburst galaxy. 



8 DISCUSSION 

Low metallicity AGN appear to be very rare objects, at least locally. 
The SDSS emission line galaxy sample only gives approximately 
40 AGN with a mass below IO'^Mq and a significant detectable 
difference in their abundances relative to typical Seyfert galaxies. 

The use of a low mass selection criterion is supported by the 
observations which reveal a correlation between mass and metal 
sensitive line ratios such as [N n]/Ha. An exact determination of a 
mass-metallicity relationship in AGN is hindered by the difficulty 
of measuring with any accuracy the metallicity in the AGN hosts. 

The 40 low metallicity AGN selected out of the 23000 Seyfert 
2s in our original SDSS sample appear to have metallicities around 
half that of typical AGN. This leads to two questions: where are the 
low metallicity AGN and why are there no AGN with even lower 
metallicities? 

We do lose some objects because we use the [O n]/13626A line 
as a diagnostic, which means we miss low metallicity objects below 
z ~ 0.03. Approximately 50 low mass AGN lie below this redshift, 
but we estimate that only 10 would fall into the low metallicity cat- 
egory. Another possibility is that low metallicity AGN are "hidden" 
by strong star formation. As discussed in[Kauffinann et al.U2003hl) . 
a large fraction of low mass galaxies are experiencing bursts of 
starformation and show strong line emission. However figure 6 in 
iKauflrmann e75] <2003d) demonstrates that it is not easy to hide a 
typical AGN of high metallicity. These authors added AGN of dif- 
ferent [O m] luminosities to low mass starburst galaxies and showed 
that at [O m] luminosities above 10'' Lq, all of them would have been 
detected. We have repeated the same analysis for low metallicity 
AGN and we find very similar results. 

It is thus clear that low metallicity AGN are rare in the 
local Universe, What about at high redshift? Our emission- 
line diagnostics involve the rest-frame optical band, and the re- 
quired spectroscopy is c hallenging at hig h redshift. Nevertheless, 
IShaplev e t al 1 2005") and Erb et al. (2006) have recently measured 
the H/?, [Om]i5007, Ha, and [Nii]i6584 emission fine fluxes in 
two small samples of galaxies at z ~ 1.4 and z ~ 2.3 respectively. 
We have plotted the locations of these galaxies in Figure^ Nearly 
all the objects lie between the Kauffmann and Kewley classification 
lines. This region is dominated by composite objects, where both 
an AGN and star formation contribute significantly to the emission- 
lines. However, the high redshift galaxies are displaced to signifi- 
cantly lower [N nj/Httratios than typical SDSS composite galaxies. 

The most natural interpretation of this is that these high red- 
shift galaxies are indeed composite galaxies, but of lower metallic- 
ity than the typical SDSS galaxies. Similar objects comprise only 
a small fraction of the SDSS sample, but in the region where most 
of the high redshift galaxies are located, bounded by the two clas- 
sification curves and the main Seyfert branch (log([Oni]/H/3) > 
3 log([Nii]/HQ')), we stifl find about 500 SDSS galaxies. Modelling 
these as composites implies that the AGN would produce 50 to 70% 
of the [Om]5007 flux, but only about 20% of the Ha and H/3 flux. 
A definitive test of the presence of an AGN would require the de- 



tection of the He n/14686 or [Ne v]/i3426 lines. Unfortunately, the 
former line is expected to be weak (a 20% AGN contribution to H/3 
implies He n/Hy8= 0.05), while the latter line is not redshifted into 
the SDSS spectral band for the great majority of objects. 

While we can not definitively say that typical SDSS galaxies 
with spectra similar to the high redshift galaxies have AGN, there 
are indirect pieces of evidence that strongly suggest that this is the 
case. When compared with a sample of nearby pure star-forming 
galaxies (all galaxies within +0.05 dex of [-0.55,0.10] on the Nii 
BPT diagram), we find that the host galaxies of these candidate 
low-metallicity composites differ in a way that is similar to what 
is found for the host galaxies of AGN ( Kauffmann et al. 2003^; 
iHeckman et"^l2004 . They have on average significantly higher 
stellar masses (~ 0.3 dex), older stellar populations (D4000 of 1.29 
compared to 1.16), redder colors (~ 0.2 greater in g - r) and weaker 
Hot equivalent widths than the pure star forming galaxies. Similar 
differences are found in their surface mass density, tr.. Concen- 
tration index and other galactic parameters. A comparison of the 
composites with objects just below the Kauffmann classification 
curve shows similar, albeit weaker, differences. However within the 
composite region there are some galaxies that are likely to be pure 
starforming galaxies as they have strong Ha equivalent widths and 
host prop erties that are sim ilar to ordinary star-forming galaxies. 

The lErbetai]<2006h object at log([Nn]/Ha) ^ -1.0 is the 
clear exception to the previous objects. Lying above both the Kauff- 
mann and Kewley classification curves, it requires a much greater 
contribution by an AGN to be offset from the star-forming branch 
due to its low metallicity.In our SDSS AGN selected sample, we 
have no low metallicity AGN or composites with which to com- 
pare this object. If we select the few SDSS objects that lie within 
the errors of the Erb object, all are clear starbursting galaxies with 
indicators such as the neon emission lines and He n/l4686A are so 
weak, if present at all, as to indicate that the contribution of an AGN 
to the spectrum is negligible. Thus this object is most likely to be a 
low metallicity starburst. 

We conclude that most of the high redshift galaxies in Fig- 
ure are most likely to be composite objects with an AGN. For- 
tunately, the AGN will have only a m odest effect on th e [N ii]/Ha 
ratio used bv lShapIev e t al. ' (2005) and lErb et alj Eo06) to measure 
the metallicity. We estimate that the AGN contribution will increase 
[N n]/Ha by only about 0. 1 to 0.2 dex, and hence cause the metallic - 
ity to be overestimated by less than . 1 dex ^ Pettini & Pagei2004) . 

If the Shaplev et al. (2005) and lErb et alJ galaxies are 

composite AGN/starforming galaxies, this may imply that low 
metallicity AGN are more common at high redshifts. We note, how- 
ever, that investigations of the metallicities of QSOs at redshifts 
greater than 1 using both broad and narrow emission lines all ap- 
pear toshow_solar to supersolar metallicities (e.g. Dietrich et alj 
l2003tlNagao et all2006l) . 

The scarcity of low metallicity AGN locally is very sugges- 
tive. The results indicate that it is unlikely for a strong AGN to be 
found in a galaxy whose mass is less than ~ IO'^Mq, or whose 
metallicity is less than solar. It is not clear whether this indicates 
that black holes only form in galaxies above some mass thresh- 
hold. There is considerable evidence that high luminosity AGN 
are hosted by gala xies with young stellar p opulations and strong 
post-burst features iKauffmann et all2003ct) . This suggests that the 
AGN must be preceded by strong star formation and this may also 
explain the high metallicities in the NLR gas. Furth er studies of 
AGN in low mass hosts (see e.g. lGreene & Hol2004h . and at early 
cosmological epochs when the mean star formation rate and gas 
densities were much higher will shed more light on these questions. 
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